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Abstract 

Theoretical foundation and application of the generalized spin-fermion 
(sp—d) exchange lattice model to magnetic and diluted magnetic semiconduc- 
tors are discussed. The capabilities of the model to describe spin quasi-particle 
spectra are investigated. The main emphasis is made on the dynamic behavior 
of two interacting subsystems, the localized spins and spin density of itiner- 
ant carriers. A nonperturbative many-body approach, the Irreducible Green 
Functions (IGF) method, is used to describe the quasi-particle dynamics. 
Scattering states are investigated and three branches of magnetic excitations 
are calculated in the regime, characteristic of a magnetic semiconductor. For 
a simplified version of the model ( Kondo lattice model ) we study the spectra 
of quasi-particle excitations with special attention given to diluted magnetic 
semiconductors. For this, to include the effects of disorder, modified mean 
fields are determined self-consistently. The role of the Coulomb correlation 

and exchange is clarified by comparing of both the cases. 
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I. INTRODUCTION 



There is considerable current interest in the properties of magnetic [1,2] and diluted 
magnetic semiconductors [3]- [8]. This field is very active and there are many aspects to the 
problem. A lot of materials were synthesized and tested [9]- [12]. The new materials design 
approach to fabrication of new functional diluted magnetic semiconductors (DMS) resulted 
in producing a variety of compounds . Diluted magnetic semiconductors are semiconducting 
alloys whose lattice contains magnetic atoms as randomly distributed substitutional impu- 
rities such as Mn-doped InAs or GaAs ( general formula A[ I J x Mn x B v ). A fraction of A 
sublattice which is substituted at random by Mn changes the carrier density from the low 
doping (nearly insulating regime) to the large doping (metallic regime ). The presence of 
the spin degree of freedom in DMS may lead to a new semiconductor spin electronics which 
will combine the advantages of the semiconducting devices with the new features due to the 
possibilities of controlling the magnetic state. 

However, the coexistence of ferromagnetism and semiconducting properties in these com- 
pounds require a suitable theoretical model which would describe well both the magnetic 
cooperative behavior and semiconducting properties as well as a rich field of interplay be- 
tween them. Existence and properties of localized and itinerant magnetism in metals, oxides 
and alloys and their interplay is an interesting but not yet fully understood problem of quan- 
tum theory of magnetism [13], [14]. Behavior and the true nature of the electronic and spin 
states, and their quasi-particle dynamics are of central importance to the understanding of 
physics of correlated systems such as magnetism and Mott-Hubbard metal-insulator transi- 
tion in metals and oxides, magnetism and heavy fermions (HF) in rare-earths compounds, 
and anomalous transport properties in perovskite manganites. This class of systems is char- 
acterized by the complex, many-branch spectra of elementary excitations. Moreover, the 
correlations effects are essential [15], [16]. 

Recently, there has been considerable interest in identifying the microscopic origin of quasi- 
particle states in such systems and a few model approaches have been proposed. Many 
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magnetic and electronic properties of magnetic semiconductors and related materials may 
be interpreted reasonably in terms of combined spin-fermion models ( SFM ) which include 
the interacting spin and charge subsystems [13], [14], [15]. This approach permits one to 
describe significant and interesting physics, e.g., the bound states and magnetic polarons 
[17], [18], anomalous transport properties, etc. 

The problem of adequate physical description within various types of generalized spin- 
fermion model has intensively been studied during the last decades, especially in the context 
of magnetic and transport properties of rare-earth and transition metals and their com- 
pounds [29], [19] and magnetic semiconductors [15]- [18], [20]- [25]. For DMS the continuous 
analog of the Zener model was proposed recently by Dietl et al. [4]. The interplay of mag- 
netic and semiconducting behavior was considered within a phenomenological mean field 
scheme. 

The majority of theoretical papers on DMS studied its properties mainly within the mean 
field approximation and continuous media terms. In such a picture the disorder effects, 
which play an essential role [30]- [32], can be taken into account roughly only. However, 
many experimental investigations call for a better understanding of the relevant physics and 
nature of solutions (especially magnetic) in terms of the lattice spin-fermion model [19], [33]. 
In the previous papers, we set up the formalism of the method of Irreducible Green Functions 
(IGF) [16]. This IGF method allows one to describe the quasi-particle spectra with damping 
for a many-particle system on a lattice with complex spectra and a strong correlation in a 
very general and natural way. This scheme differs from the traditional method of decoupling 
of an infinite chain of equations [26] and permits construction of the relevant dynamical so- 
lutions in a self-consistent way at the level of the Dyson equation without decoupling the 
chain of equations of motion for the GFs. 

In this paper, we apply the IGF formalism to consider the quasi-particle spectra for the 
lattice spin-fermion model consisting of two interacting subsystems. It is the purpose of this 
paper to explore more fully the notion of Generalized Mean Fields (GMF) [16] which may 
arise in the system of interacting localized spins ( including effects of disorder ) and lattice 
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fermions to justify and understand the "nature" of the relevant mean fields. Theoretical 
foundation and application of the generalized spin-fermion ( sp-d ) exchange model to mag- 
netic and diluted magnetic semiconductors are discussed in some detail. The capabilities of 
the model to describe quasi-particle spectra are investigated. The purpose is to investigate 
the quasi-particle spectra and GMF of the magnetic semiconductors consisting of two in- 
teracting charge and spin subsystems within the lattice spin-fermion model in a unified and 
coherent fashion to analyze the role and influence of the Coulomb correlation and exchange. 



II. SPIN-FERMION MODEL 

The concept of the sp — d { or d — f ) model plays an important role in the quantum 
theory of magnetism [13], [14], [19], [27], [28], [29]. In this section, we consider the generalized 
sp — d model which describes the localized 3d(4/)-spins interacting with s(p)-like conduction 
(itinerant) electrons ( or holes ) and takes into consideration the electron-electron interaction. 
The total Hamiltonian of the model is given by 

H = H S + H s _ d + H d (1) 

The Hamiltonian of band electrons (or holes) is given by 

H » = 12 ^ tiAo-i" + \ 12 U n * n i-° ( 2 ) 
ij <y io 

This is the Hubbard model. We adopt the notation 

a ia = N~ 1/2 a ka exp(ikRi) a\ a = N~ 1/2 ^ a\ a exp(-ikRi) 
k k 

In the case of a pure semiconductor, at low temperatures the conduction electron band is 
empty and the Coulomb term U is therefore not so important. A partial occupation of the 
band leads to an increase in the role of the Coulomb correlation. It is clear that we treat 
conduction electrons as s-electrons in the Wannier representation. In doped DMS the carrier 
system is the valence band p-holes. 
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The band energy of Bloch electrons e(k) is denned as follows: 

*« = iV_1 E exp[ik(Ri - Rj)}, 
k 

where N is the number of lattice sites. For the tight-binding electrons in a cubic lattice we 
use the standard expression for the dispersion 

e(k) = 2^2 *(0 cos(A;a a ) (3) 

a 

, where a a denotes the lattice vectors in a simple lattice with the inversion centre. 

The term H s ^d describes the interaction of the total 3d(4f)-spin with the spin density of the 

itinerant carriers [28] 

H s _ d = -2 £ IBiSi = -IN- 1 / 2 EE^oLo^ + z a Si q ala k+qa ] (4) 

i kq a 

where sign factor z c is given by 

z a = (+or— ) /or o- = (| or |) 

and 

(SZ Q ifa = + 

sz: = I 

S+ if a = - 



In DMS the local exchange coupling resulted from the p — d hybridization between the Mn 
d levels and the p valence band / ~ V p 2 _ d ■ For the subsystem of localized spins we have 

Hd = —77 E JijSiSj = ~aE JqS q S- q (5) 
Z ij Z 1 

Here we use the notation 

S? = N' 1 ' 2 £ St expiikRi) St = N- 1 ' 2 ]T St exp(-^) 

2 

[Sk > ^q] = T'S'fe+g , ^ ] = j^i/ 2 ^k+q 

^3 = N' 1 E ^jfcexp[ifc(fli - 



This term describes the direct exchange interaction between the localized 3d (4f) magnetic 
moments at the lattice sites i and j. In the DMS system this interaction is rather small. The 
ferromagnetic interaction between the local Mn moments is mediated by the real itinerant 
carriers in the valence band of the host semiconductor material. The carrier polarization 
produces the RKKY exchange interaction of Mn local moments 

Hrkky = — £ KijSiSj (6) 

We emphasize that ~ \P\ ~ V^_ d . To explain this, let us remind that the microscopic 
model [28], which contains the basic physics, is the Anderson-Kondo model 

# = £ £ tiAaaj* ~ V J2 J2( a ta d j« + h - c -) 

ij a ij a 

-^££^t + ^£^UU (7) 

i cr i a 

For the symmetric case U = 1E& and for U ^> V Eq.(7) can be mapped onto the Kondo 
lattice model ( KLM ) 

H = £ £ t ija la ja - £ (8) 

ij a i 

Here I ~ %■ The KLM may be viewed as the low-energy sector of the initial model Eq.( 
7). 

III. OUTLINE OF THE IGF METHOD 

In this section, we discuss the main ideas of the IGF approach that allows one to describe 
completely quasi-particle spectra with damping in a very natural way. 
We reformulated the two-time GF method [16] to the form which is especially adjusted to 
correlated fermion systems on a lattice and systems with complex spectra [15]- [19]. A very 
important concept of the whole method is the Generalized Mean Fields (GMFs), as it was 
formulated in [16]. These GMFs have a complicated structure for the strongly correlated 
case and complex spectra, and are not reduced to the functional of mean densities of the 
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electrons or spins when one calculates excitation spectra at finite temperatures. 
To clarify the foregoing, let us consider a retarded GF of the form 



G r =« A(t), A\t') »= -i6(t - t') < [A(t)A\t% >, 7] = ±1 (9) 

As an introduction to the concept of IGFs, let us describe the main ideas of this approach 
in a symbolic and simplified form. To calculate the retarded GF G(t — t'), let us write down 
the equation of motion for it 

uG(lo) =< [A, A% > + « [A, H\_ | A ] » UJ (10) 

The essence of the method is as follows [16]: 

It is based on the notion of the "IRREDUCIBLE" parts of GFs (or the irreducible parts of 
the operators, A and A\ out of which the GF is constructed) in terms of which it is possible, 
without recourse to a truncation of the hierarchy of equations for the GFs, to write down 
the exact Dyson equation and to obtain an exact analytic representation for the self-energy 
operator. By definition, we introduce the irreducible part (ir) of the GF 

(ir) « [A,H\_\rf »=« [A,H}_ - zA\ti » (11) 

The unknown constant z is defined by the condition (or constraint) 

<[[A,Ht\AX>=0 (12) 

which is an analogue of the orthogonality condition in the Mori formalism. From the condi- 
tion (12) one can find 

<[[A,H].,A% > Mi 

<[A,A%> M [i6) 

Here M and Mi are the zeroth and first order moments of the spectral density. Therefore, 
the irreducible GFs are defined so that they cannot be reduced to the lower-order ones by 
any kind of decoupling. It is worth noting that the term "irreducible" in a group theory 
means a representation of a symmetry operation that cannot be expressed in terms of lower 



dimensional representations. Irreducible (or connected ) correlation functions are known in 
statistical mechanics. In the diagrammatic approach, the irreducible vertices are defined 
as graphs that do not contain inner parts connected by the G°-line. With the aid of the 
definition (11) these concepts are expressed in terms of retarded and advanced GFs. The 
procedure extracts all relevant (for the problem under consideration) mean-field contribu- 
tions and puts them into the generalized mean-field GF which is defined here as 

^ = <[A,A%> 
(u - z) 

To calculate the IGF (ir) << [A, H]_(t), A^(t') » in (10), we have to write the equation of 
motion for it after differentiation with respect to the second time variable t'. The condition 
of orthogonality (12) removes the inhomogeneous term from this equation and is a very 
crucial point of the whole approach. If one introduces the irreducible part for the right-hand 
side operator as discussed above for the "left" operator, the equation of motion (10) can be 
exactly rewritten in the following form 

G = G° + G°PG° (15) 

The scattering operator P is given by 

P=(M )- 1 ( (ir) « [A,HU[A\H}_ »^)(M )- 1 (16) 

The structure of equation ( 16) enables us to determine the self-energy operator M, by 
analogy with the diagram technique 

P = M + MG°P (17) 

From the definition (17) it follows that the self-energy operator M is defined as a proper 
(in the diagrammatic language, "connected") part of the scattering operator M = (P) p . 
As a result, we obtain the exact Dyson equation for the thermodynamic double-time Green 
functions 

G = G° + G°MG (18) 



The difference between P and M can be regarded as two different solutions of two integral 
equations (15) and (18). But from the Dyson equation (18) only the full GF is seen to be 
expressed as a formal solution of the form 

G = [(G )- 1 - M}- 1 (19) 

Equation (19) can be regarded as an alternative form of the Dyson equation (18) and the 
definition of M provided that the generalized mean-field GF G° is specified. On the contrary 
, for the scattering operator P, instead of the property G°G~ 1 + G°M = 1, one has the 
property 

(G )- 1 - G- 1 = PG°G- 1 

Thus, the very functional form of the formal solution (19) determines the difference between 
P and M precisely. 

Thus, by introducing irreducible parts of GF (or irreducible parts of the operators, out 
of which the GF is constructed) the equation of motion (10) for the GF can exactly be ( 
but using the orthogonality constraint (12)) transformed into the Dyson equation for the 
double-time thermal GF (18). This result is very remarkable because the traditional form 
of the GF method does not include this point. Notice that all quantities thus considered 
are exact. Approximations can be generated not by truncating the set of coupled equations 
of motions but by a specific approximation of the functional form of the mass operator M 
within a self-consistent scheme expressing M in terms of initial GF 

M « F[G] 

Different approximations are relevant to different physical situations. 

The projection operator technique has essentially the same philosophy. But with using the 
constraint (12) in our approach we emphasize the fundamental and central role of the Dyson 
equation for the calculation of single-particle properties of many-body systems. The problem 
of reducing the whole hierarchy of equations involving higher-order GFs by a coupled nonlin- 
ear set of integro-differential equations connecting the single-particle GF to the self-energy 
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operator is rather nontrivial. A characteristic feature of these equations is that, besides the 
single-particle GF, they involve also higher-order GF. The irreducible counterparts of the 
GFs, vertex functions, serve to identify correctly the self-energy as 

M = Gq 1 - G- 1 

The integral form of Dyson equation (18) gives M the physical meaning of a nonlocal and 
energy-dependent effective single-particle potential. This meaning can be verified for the 
exact self-energy through the diagrammatic expansion for the causal GF. 
It is important to note that for the retarded and advanced GFs, the notion of the proper 
part M = (P) p is symbolic in nature [16]. In a certain sense, it is possible to say that it 
is defined here by analogy with the irreducible many-particle T-matrix. Furthermore, by 
analogy with the diagrammatic technique, we can also introduce the proper part defined 
as a solution to the integral equation (17). These analogues allow us to understand better 
the formal structure of the Dyson equation for the double-time thermal GF but only in a 
symbolic form . However, because of the identical form of the equations for GFs for all three 
types ( advanced, retarded, and causal ), we can convert in each stage of calculations to 
causal GF and, thereby, confirm the substantiated nature of definition (17)! We therefore 
should speak of an analogy of the Dyson equation. Hereafter, we drop this stipulating, since 
it does not cause any misunderstanding. In a sense, the IGF method is a variant of the 
Gram-Schmidt orthogonalization procedure. 

It should be emphasized that the scheme presented above gives just a general idea of the IGF 
method. A more exact explanation why one should not introduce the approximation already 
in P, instead of having to work out M, is given below when working out the application of 
the method to specific problems. 

The general philosophy of the IGF method is in the separation and identification of elastic 
scattering effects and inelastic ones. This latter point is quite often underestimated, and 
both effects are mixed. However, as far as the right definition of quasi-particle damping is 
concerned, the separation of elastic and inelastic scattering processes is believed to be cru- 
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daily important for many-body systems with complicated spectra and strong interaction. 
From a technical point of view, the elastic GMF renormalizations can exhibit quite a non- 
trivial structure. To obtain this structure correctly, one should construct the full GF from 
the complete algebra of relevant operators and develop a special projection procedure for 
higher-order GFs, in accordance with a given algebra. Then a natural question arises how 
to select the relevant set of operators {Ai, A 2 , ...A n } describing the "relevant degrees of 
freedom". The above consideration suggests an intuitive and heuristic way to the suitable 
procedure as arising from an infinite chain of equations of motion (10). Let us consider the 
column 

(AA 

A 2 



where 



Al =A, A 2 = [A, H] , A 3 = [[A,H],H],...A n = [[...[A, H]...H] 



Then the most general possible Green function can be expressed as a matrix 

A 2 



G=« 



. Ai A\ 



At) » 



\A n ) 



This generalized Green function describes the one-, two- and n-particle dynamics. The equa- 
tion of motion for it includes, as a particular case, the Dyson equation for single-particle 
Green function, and the Bethe-Salpeter equation which is the equation of motion for the 
two-particle Green function and which is an analogue of the Dyson equation, etc . The 
corresponding reduced equations should be extracted from the equation of motion for the 
generalized GF with the aid of the special techniques such as the projection method and 
similar techniques. This must be a final goal towards a real understanding of the true many- 
body dynamics. At this point, it is worthwhile to underline that the above discussion is a 
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heuristic scheme only but not a straightforward recipe. The specific method of introducing 
the IGFs depends on the form of operators A n , the type of the Hamiltonian, and conditions 
of the problem. 

Here a sketchy form of the IGF method is presented. The aim is to introduce the general 
scheme and to lay the groundwork for generalizations. We demonstrated in [16] that the 
IGF method is a powerful tool for describing the quasi-particle excitation spectra, allowing 
a deeper understanding of elastic and inelastic quasi-particle scattering effects and the cor- 
responding aspects of damping and finite lifetimes. In the present context, it provides an 
efficient tool for analysis of the mean fields and generalized mean fields of the complicated 
many-body models. 

IV. QUASI-PARTICLE DYNAMICS OF THE (SP - D) MODEL 

To describe self-consistently the spin dynamics of the extended sp — d model, one should 
take into account the full algebra of relevant operators of the suitable "spin modes" which 
are appropriate when the goal is to describe self-consistently the quasi-particle spectra of 
two interacting subsystem. 
We have two kinds of spin variables 

Ski S-k = i^kV 

a k =Y1 a tf a k+qh a -k = ( a k ) f = H a k+ql a q1 

q q 
Let us consider the equations of motion 

[Ski Hs-d]- = -IN 1 ^2[2S^_ q a p ^a p+q i - S^_ q (a p ^a p+q i - a pl a p+q i)] (20) 

pq 

[S_ k , H s - d ]- = -IN' 1 X][ 2 ^fc-g a pJ. a P+9T ~ ^fc- g ( a pT a P+9-T ~ a ll a p+ql)) ( 21 ) 

pq 

[s z k , = -in- 1 Y,(SLAl a p + q1 ~ s*"- 9 4tw) ( 22 ) 



pq 
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[S+, H d ]_ = N-^ £ US z q St q - S z k _ q S+) (23) 

[SZ k , H d \_ = N- 1 ' 2 £ J q (S z - {k+q) S; - S*SZ {k+q) ) (24) 
[4t a <?+fc|> H s ]_ = (e(q + k) — e(q))a^a q+kl + 

U N' 1 Y( a l-l a l+p'^ a p\ a q+p'+kl - a l+p'^ a l-p'i a pl a q+k.l) ( 2 5) 



PP 



[al^a q+kl , # s -<i]- - IN 1/2 Y^S-p'( a i\ a p+p'^p,<i+k ~ a li a q+kiS q ,p+ P ') 

-Sl pl (a q ^a p+P > iS Pt<1+ k + a p ^a q+ ki5 qt p +p ')] (26) 



pp 



From Eq.(20) - Eq.(26) it follows that the localized and itinerant spin variables are coupled. 
Suitable algebra of relevant operators should be described by the 'spinor' (j^J ("relevant 
degrees of freedom"), according to the IGF strategy. 

The model Hamiltonian H = H s + H s _ d + Hd was used in [34], [35] for calculations of the 
spin-wave spectra and was called the modified Zener model. In this model, as applied to 
transition metals, the itinerant electrons are described by a Hubbard Hamiltonian and the 
itinerant electron couples the localized spin (Hund's rule coupling) by a term H s _d- Because 
of the inequivalent spin systems, localized and itinerant, a consequence of the model is the 
existence of acoustic and optic branches of the quasi-particle spectrum of spin excitations. In 
DMS the local antiferromagnetic interaction H s _d produces the coupling between the carriers 
( which are holes in GaMnAs ) and the Mn magnetic moments ( s = 5/2) which leads to 
ferromagnetic ordering of Mn spins in a certain range of concentration. The Kondo physics 
is irrelevant in this case, but the fully determined and consistent microscopic mechanism of 
the ferromagnetic ordering is still under debates [3]. An important question in this context 
is the self-consistent picture of the quasi-particle many-body dynamics which takes into 
account the complex structure of the spectra. 
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A. Spin Dynamics of the s — d Model. Scattering Regime. 



In this section, we discuss the spectrum of spin excitations in the sp — d model. We 
consider the double-time thermal GF of localized spins [26] which is defined as 

Gr+-( fc; t - if) =« S+(t), Sz k (f) »= -i6(t - if) < [S+(t), SZ k (t')}_ >= 

/+oo 
dujexp{-iut)G + -{k;uj) (27) 
-oo 

The next step is to write down the equation of motion for the GF. Our attention will be 
focused on spin dynamics of the model. To describe self-consistent ly of the spin dynamics of 
the sp — d model, one should take into account the full algebra of relevant operators of the 
suitable " spin modes" which are appropriate when the goal is to describe self-consistently the 
quasi-particle spectra of two interacting subsystems. We introduce the generalized matrix 
GF of the form 

«S k \SZ k » «S k \aZ k »' 
« o k \SZ k » « o k \oZ k » 

Here 



G{k-u) (28) 



G'(k;u), (29) 



i q 
Equivalently, we can do the calculations with the matrix of the form 

« S k\ S -k » « S k\ a k+ql a q1 » 

« al ] a q+kl \S_ k » « af^aq+ki^+^a^ » 
but the form of Eq.(28) is slightly more convenient. 

Let us consider the equation of motion for the GF G(k;uj). By differentiation of the GF 
<< S k (t)\B(t') » with respect to the first time, t, we find 

„ « S t \B »„= { 2JV - 1/2 < S » >} + (30) 

— « S k-q( a ll a p+q]- - a tl a p+ql) ~ 2S k-q a U a P+ll\ B >>^ 



pq 



+7V-V2 ^ Jq « {S * S +_ q _ Sl_ q S+)\B » w (31) 
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where 

B 



Let us introduce by definition irreducible (ir) operators as 

(S*) %T = S q — < Sq > 5g : o; (o,p +qa a po -) %T = Q>p+ qcr Q>pa— < a p a a P<r > ^<3,0 (32) 

((S z q ) ir St q - (S' k . q ) ir Sf) fr = ((S z ) ir S+_ q ~ (S z k _ q ) *S+) - (0, - <t> k - q )St (33) 
From the condition (12) 

< [((S z q ) ir S+_ g - (S z _ q ) ir S+ - (0, - <f> k - q )s+), sz k }- >= o 

one can find 

2K z a z + K~+ , 
<P q = — „ 34 

K? =< (S z ) ir (S* q ) ir >; K q + =< SZ q S+ > (35) 

Using the definition of the irreducible parts the equation of motion Eq.(30 ) can be exactly 
transformed to the following form: 



•ATl/2 



fix « S+\B » u +Q 2 « a+\ B » u = I v 1 \+ « A^B » u (36) 



where 



<r c ! z >> 2K ZZ 4- K~ + 

"i = <" - J » - ■« - JV_1/2 - ■*--*> 2 <gg> " /( "' " " l) (37) 

V, = (38, 

1 1 

n <^ = at < a \o a i° >= ^J2 for = H( ex P(/5 e (^)) + 1) 

q q q 

e(qa) = e(q) - zJN' 1 ^ < S z > +Un^ a 

^ = E(™t + n i); 0<n<2 

The many-particle operator Ai reads 
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^! ~~ JTf^2[Sk-q( a i\ a p+qT ~~ a p| a p+?l) ~ ^i^k-q) ^ a p1 a p+ql\ 
pq 

+7V-V2 £ J^) ^ _ {SLq) ir S ^ ir ( g g) 
<? 

and it satisfies the conditions 

< [A u SZ q ]- >=< [A u aZ q ]- >=0 

To write down the equation of motion for the Fourier transform of the GF 
<< a£(t), B(t') >>, we need the following auxiliary equation of motion: 

(w + e(p) -e(p + k) -2IN~ 1/2 < S% > -C/(ra T - raj) « a^a p+kl \B » u + (40) 
UN~\f lA - f p+H ) « a+\B » u +IN- 1 / 2 (f pl - f p+kl ) « S+\B » u = 

\fr , A - IN~ l/2 « ^-r( a p-[ a q+r1^p+k,q ~ a q \ a p+k^p,q+r) W \B »u> ~ 

[Up] - Jp+kl) J qr 

IN- 1/2 Y,« (S z - r ) ir (a^a p+kl 5 p>g+r + a^a q+ri 8 p+k , q )\B » u + 

qr 

UN' 1 « (a\ }] a q+r] a q ^a p+r+kl - a ] p+r] a q _ ri a q[ a p+k {) ir \B » w 

qr 

Let us use the following notation: 



A 2 — -IN 1/2 ^2[SZ r (a P ^a q+r -i5 p+k>q - a qi a p+k i8 Pjq+r ) " 

qr 

' l q^ u P+kl u p,q+r T 

qr 

u p , k = (lu + e(p) - e(p + fc) - A) (42) 



(S! r ) ir (a q ^a p+ki 5 p>q+r + a pT a 9+ri 5 p+fc;? )] + (41) 
UN' 1 ^2{a p ^a q+r] a q ^a p+r+H - a p+r ^a q _ ri a qi a p+ki ) ir 



A = 2IN~ 1/2 <S*> -C/(ra T - raj = 21 S - Um = A/ + Ay (43) 

X -(fe,a;)^7V-^ ^|-/pT) (44) 

Now we consider the GF << <j£(t), B(t') » . Similarly to Eq.( 36), we have 

-N^Ix^Ku) « S+\B » u +(1 - Ux^k,u)) « <j+\B » u = 

I M L n}+E— «^|S» W (45) 

Here the following definition of the irreducible part for the Coulomb correlation term was 
used 
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( a ll a l+r^ a <lt a P+r+kl — a l + r]- a l-rl a (ll a P+kl) " 

= (a P ^a q+r ^a q ^a p+r+k i — a p+r ^a q _ r ^a q ia p+k i) 

t At 

- < a q _ ri a qi > 5 q - ryq a} p+r] a p+ki (46) 
The operator A 2 satisfies the conditions 

< [A 2 ,SZ k \- >=< [A 2 ,vZ k }->=0 

In the matrix notation the full equation of motion for the GF G(k; u>) can now be summarized 
in the following form: 

ClG(k;u) = i + Y,*(p)D(p;u) (47) 
p 

(G(^))t(n)t = (/)t + ^0( p;w ))t($ (p ))t 



where 



fii tt 2 \ (I- l N^ 2 n 2 

n=\ / = " I (48) 

-IN 1 ^ (1-U X S )J V -N X S . 

i\S k » <<Ai|o-Z fe >>\ ( N' 1 

D(p; u) — I S>(p) = i I (49) 



To calculate the higher order GFs in ( 47), we differentiate its r.h.s. with respect to the 
second-time variable (f). Let us give explicitly one of the four equations. After introducing 
the irreducible parts as discussed above we get 

<< A\SZ k »u> = « A \^-(k-q')( a i'f a p'+g"t - ap'l'V+9'i) ir - 25'i (A ._ ? , ) a^ i a p / +g / T >> w 

p'q' 



+n- II2 Y.J^ « A \[{s z q >) ir sz {k+ql) - (s z _ {k+q/) ) ir s q ,\ - » 

q' 



Here, the quantity A in the l.h.s. should be substituted by 

f(Qz\ ir c+ ( cz \ ir c+\ ir 

IWeJ D fe-(3 ~~ Wk-q) D o ) 



A = I 



Sk-q( a il a p+q1 ~ a ll a p+ql) %T 
2S k -q a il a p+qi 
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In the matrix notation the full equation of motion for the GF D(k;u) can now be written 
in the following form: 

nD(p;uj)=J2Hp')D 1 (p'^) (51) 
p> 

where 

(«A 1 \A\» «a 1 \a\»\ 

Di=[ , . (52) 

\«A 2 \A\» «A 2 \A\») 

Combining both (the first- and second-time differentiated) equations of motion, we get the 
"exact" ( no approximation has been made till now) "scattering" equation 

QG(k;u) = I + Y,Hp)P(p,p')Hp')(&)- 1 (53) 
pp' 

This equation can be identically transformed to the standard form Eq.( 16) 

G = G + Go^I-'HpMPiP'Mp')!- 1 )^ 

pp' 

G = G + GoPGo (54) 

Here we have introduced the generalized mean-field (GMF) GF G , according to the following 
definition: 

Go = (55) 
The scattering operator P has the form 

p = ^ l ^l»(p)p(p,p / )^(p , )^ 1 (56) 
pp' 

Here we have used the obvious notation 

(«Ai\A\» «Ai\A\»\ 
P(p,j/;u;)= (57) 
\«A 2 \A\» «A 2 \Al»J 

As is shown above, Eq.( 54) can be transformed exactly into the Dyson equation ( 18) 

G = G + G MGo (58) 
with the self-energy operator M given as 

M = (P) p (59) 
Hence, the determination of the full GF G has been reduced to that of Go and M. 
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V. GENERALIZED MEAN-FIELD GF 



From the definition ( 55), the GF matrix in the generalized mean-field approximation 
reads 

((l-uxDi-Wfish n 2 N X s \ 

G = iT 1 (60) 
V V 2 N X S -^Nx'oJ 

where 

R = (l- UxD^li + ^IN^xl 



Let us write down explicitly the diagonal matrix elements Gq 1 and Gq 2 

2S 

« SJ|S: ' >>0= a l + 2PS M (61) 

« a k \a_ k » = ni + 2pSxS{ki0) (62) 

where 

X s (k,u) = Xo(k,uj)(l - UxZfcoj))- 1 (63) 
S = N- 1/2 < S* > 

To clarify the functional structure of the generalized mean-field GFs (61 ) and (62 ), let us 
consider a few limiting cases. 

A. Uncoupled Subsystems 

To clarify the calculation of quasiparticle spectra of coupled localized and itinerant sub- 
systems, it is instructive to consider an artificial limit of uncoupled subsystems. We then 
assume that the local exchange parameter 1 = 0. In this limiting case we have 

2S 

« S+\S_ k »°= - _ F{ Jq _ Jfe) _ _i_ Eq(Jq _ Jq _ k) (2K? + K-+) (64) 
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«a+\aI k »°= X s (k,u;) (65) 

The spectrum of quasi-particle excitations of localized spins without damping follows from 
the poles of the generalized mean-field GF (64) 

u(k) = S(J ~Jk) + ~ £(■/, - J q -k)(2K» + K~+) (66) 

It is seen that due to the correct definition of generalized mean fields we get the result for 
the localized spin Heisenberg subsystem which includes both the simplest spin-wave result 
and the result of Tyablikov decoupling as limiting cases. In the hydrodynamic limit k — > 0, 
uo — > it leads to the dispersion law u(k) = Dk 2 . 
The exchange integral Jk can be written in the following way: 

J fc = 5>xp(-M,)J(|^|) (67) 

i 

The expansion in small k gives 

J k = £ J(|^|) - ^E(^)V(|^|) = Jo -|E(*(I^I) (68) 

i i i 

Here n = k/k is the unit vector. The values Jh- q can be evaluated in a similar way 

J k _ q — Jq— (kV q )J q + ^(kW q ) 2 J q + ■■■ (69) 

(kV q )J q = -iJ20Ri)J(\Ri\)exp(-tqR t ) 

i 

{kVqfJq = -W(^) 2 A\Ri\)^P(-iqRi) 

i 

Combining Eq.(69), Eq.(68), and Eq.(66) we get 

2 9 

« St\SZ k »°= 7T- (70) 

fel ~ k u-uj(k) 

u{k - 0) = (S(J ~Jk) + ~ £(■/, - J 9 -*)(2^* + K: + )) ~ -Di/c 2 



2A^S* 

1 2 ' " 2,S'- ? 



C AT 

9 
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Let us now consider the spin susceptibility of itinerant carriers Eq. (65 ) in the hydrodynamic 
limit k — > 0, uj — > 0. It is convenient to consider the static limit of Eq. (65 ) 



«a k \a_ k » U= 1 _ t/xS(fcj0) (71) 

x5(*> o) = 1 E <q) ^ + ^_ Au 



= U(ni — ri|) = Um 
To proceed, we make a small-k expansion of the form 



e ( 5 + fc) - e(q) = (kV q )e(q) + ^(W g ) 2 e(g) + • • • (72) 



N A*Y J " J *»2" ™ NA* q 
The poles of the spin susceptibility of itinerant carriers are determined by the equation 

1-Ux s (k,cu) = (73) 

In another form this reads in detail 

1 = ^_ sr (si ~ fq+H 

! " N Y e(k + q) - e(q) + A v - uj 

If we set uj = E{k) and put then k = 0, we get the equation for the excitation energy 
E(k = 0) 



N ^ Au - E(k = 0) -£(£; = 0) C/ 

which is satisfied if E(k = 0) = 0. Thus, a solution of Eq.(73 ) exists which has the property 
limfc^o E{k) = and this solution corresponds to an acoustic spin-wave branch of excitations 

E(k) = D 2 k 2 = -^-;E(/«T + Ul)(kV q )h(q) + ^|rE(/,T - U)(^AQ)) 2 (74) 

u = e(k + q) - e(q) + A v 

It is seen that the stiffness constant D% can be interpreted as expanded in For the 
tight-binding electrons in s.c. lattice the spin wave dispersion relation D 2 k 2 becomes 
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D 2 e = (S^-n^B ^V ^W ^T - ^y^ V;e(g)] = (75) 
2t 2 a 2 

(3(ra T - nj) — — 5Z(/«T _ sin(g x a) + fc y sin(g y a) + k z sm(q z a)) 2 - 

ta2 ^2U'q1 + fql)(kl cos q x a + k 2 cos q y a + k\ cos q z a)) 



B. Coupled Subsystems 

The next stage in the analysis of the quasi-particle spectra of the (sp — d) model is the 
introduction of the nonzero coupling I. The full generalized mean field GFs can be rewritten 

as 

2S 



« StlS'u »°= 



uj-Im- S(J -J k )-zks - J i-^ K l Z + K ^ + ) + 2p Sx s (k, u) 

(76) 



« a+\a- h »°= - t^t^TTi ^ (77) 
Here the notation is used 

2I 2 S 

Ueff = U ~ Lu-Im ] m= ^- n i) 

The expression Eq.(77) coincides with that for the itinerant spin susceptibility as calculated 
in [13]. It is instructive to consider separately the four different cases, 

(i) / ^0,J = 0,[/ = 

(ii) f^o,J^O,£/ = 

(iii) I^0,J = 0,U^0 

(iv) J^0,J^0,£/^0 
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The first case I ^ 0, J = 0,U = corresponds to a model which is commonly called the 
Kondo lattice model. It can be seen that GFs (76 ) and (77 ) are then equal to 

2 9 

« S+\SZ k »°= = or9 ^ ~rj — r (78) 



« <j+|<7l fc » u = (79) 

In order to calculate the acoustic pole of the GF (78 ), we make use of the small (k,u) 
expansion. Hence we get 

« S+\SZ k (80) 

2g(l + 1)- 1 

w " (! + D" 1 ^ + /,i)(^V ff )M9) - lAj E,(/, T - / ff i)(^V ff e(9))2] 

It follows from Eq.(80 ) that the stiffness constant D is proportional to the total magneti- 
zation of the system. 



2. 



In the second case I ^ 0, J ^ 0,U = 0, we get 

«s+\s-_ k »°- 



uo-Im- S(J -J k )-^ks £<A " J q -k)(2K% z + K-+) + 2PS X s (k, u) 

(81) 



« at\aZ k » u = ZY s(h \ ( 82 ) 

In order to calculate the acoustic pole of the GF (81 ), we make use of the small (k,u) 
expansion again. We then get 



« s+\sz k » ^ 



2^1 + j)- 1 



u - (1 + i)- 1 ^ 2 - (1 + f)- 1 ^ + f q i)(kV q ye(q) - ^E ? (/ 9T - /,i)(*V ff e(9)) s 

It follows from Eqs.(80 ) and (83 ) that the stiffness constant D is proportional to the total 
magnetization of the system. 
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3. 



The third case I ^ 0, J = 0,U ^ corresponds to a model which is called the modified 
Zener lattice model [34]. It can be seen that in this case GFs (76 ) and (77 ) are equal to 

2 9 

« S+\SZ k »°= ; or9 ^ „ r (84) 

« *t\a- k »°= - T^r^-n ; (85) 

The results obtained here coincide with those of Bartel [34]. The excitation energies for 
the localized spin and spin densities of itinerant carriers are found from the zeros of the 
denominators of << S£\SZ k »° and << \aZ k »° which yield identical excitation 
spectra, consisting of three branches, the acoustic spin wave E ac (k), the optical spin wave 
E op (k), and the Stoner continuum E st (k) 

E ac (k) = Dk 2 

E°i>(k) = E*>-D(l-—)k 2 ; E° p = I(m + 2S) 
E st (k) = e(k + q) — e(q) + A 



I 

The most general is the forth case, I ^ 0, J ^ 0, U ^ 0. The total GF of the coupled 
system is given by Eq.(76 ). The magnetic excitation spectrum follows from the poles of the 
GF ( 60) 

R = (i-u x s )n 1 + n 2 iN 1 / 2 x s = o 

and consists of three branches - the acoustic spin wave E ac (k), the optical spin wave E op (k), 
and the Stoner continuum E st (k). 

Let us, as a first approximation, consider the last term in its denominator , which is the 
dynamic spin susceptibility of itinerant carriers, in the static limit, without any frequency 
dependence. The GF Eq.(76 ) then becomes equal to 
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C+ 1 c~ o~ 

« k\ -fc >> ~ - _ Im _ s{Jo _ Jk) _ _i_ ^-p- _ Jq _ k){ 2 K ,, + + 2 J2^(A;, 0) 

(86) 

It is possible to verify that in the limit /c — > 

2/ 2 % s (A;,0) » Im - ^(/, t + / gi )(W g ) 2 e(g) + E(/«T " Al)(^e(g)) 2 (87) 

Then for a;, — >■ Eq.(87 ) becomes 

« St\SZ k »°^ (88) 

2S 

oj - Dtf - + /,i)(fcV 9 ) 2 e(g) + £ 9 (/*T ~ /,i)(*V,e(g))2 

This expression can be expected to be qualitatively correct in spite of the primitive approx- 
imation. The spectrum of Stoner excitations is given by 

E s \k) = e(k + q) - e(q) + A (89) 

In addition to the acoustic branch there is an optical branch of spin excitations. This can 
be seen from the following: For k = we get for R = the quadratic equation in uj with 
two solutions, uj = and uj = I(m + 2S) = Eq P . In the hydrodynamic limit, k — > 0, uj — > 
the GF (75) can be written as 

2S 

« S£\SZ k » ^ f^-ttt (90) 

k 1 k uj- E ac (k) v 1 

where the acoustic spin wave energies are given by 

E ac {k) = Dk 2 = (|[^o + iEW + (91) 



, 2 Lru 2Ar52 
1 1 „ . , _ 1 1 



E(AT + /,i)(*V ff ) 2 *(9) + tt^EC/W - / 9 |)(nV g e(g)) 2 )A; 2 



2A^2^Y iVA2S „ 

For the optical spin wave branch the estimations can be carried out as in paper [13] 

E op (k) = E° P - D op k 2 (92) 

In the GMF approximation the density of itinerant electrons ( and the band splitting A) 
can be evaluated by solving the equation 
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n ° = -Jt £[exp(/3(e(£;) + Un_ a -IS-e F )) + l]' 1 (93) 

k 

Hence, the stiffness constant D can be expressed by the parameters of the sp — d model 
Hamiltonian . 

VI. EFFECTS OF DISORDER IN DMS 

We now proceed to give a simple and qualitative discussion of the effects of disorder in 
DMS. The main aim of the investigation of DMS is to give a successful microscopic picture of 
the ferromagnetic ordering of localized spins induced by the interaction with the spin density 
of itinerant charge carriers. The big amount of information was accumulated on the problem 
of ordering of localized spins of Mn 2+ ions in Cdi- x Mn x Te and similar compounds [1], [37]. 
Recently, main attention was directed to studying III-V semiconducting alloys doped with 
Mn. As it was stated above, a suitable model, which may be used for investigation of this 
problem ( at least at the initial stage ) is a modified Kondo lattice model (8 ) 

H = Y1J2 UjoLaja - J2 2Iiy i°iSi ( 94 ) 

ij o i 

Here ui projects out sites occupied by Mn atoms, i.e.: 




1 if site i is occupied by Mn 
if site i is occupied by Ga 



This model is relevant for the doped II- VI or III-V compound. The essential feature of the 
model is that it describes a mechanism of how the spins of carriers (electrons or holes ) be- 
come polarized due to the local antiferromagnetic exchange interactions with localized spins. 
In A[ I J x Mn x B v the main magnetic interaction is an antiferromagnetic exchange between the 

-, — * — * 

Mn spins and the charge-carrier spins. The superexchange term = JijSiSj is an- 

tiferromagnetic also but is as a rule rather small in the concentration range of interest ( 
x ~ 0.05). In the case of Mn-doped III-V compounds the antiferromagnetic superexchange 
interaction will generally reduce the ferromagnetic ordering temperature. As a result, the 
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carrier-induced ferromagnetism in DMS may arise due to the effective ferromagnetic inter- 
action between the Mn spins. In other words, the ferromagnetism in this system is most 
probably related to the uncompensated Mn spins and is mediated by holes. The density 
of Mn ions cun is greater than the hole density p, CMn P- The optimal interrelation of 
both the magnitudes is a delicate and subtle question and was analyzed recently in paper 
[36]. It was shown that the concentration of free holes and ferromagnetically active Mn 
spins were governed by the position of the Fermi level which controls the formation energy 
of compensating interstitial Mn donors. The experimental evidence has been provided that 
the upper limit of the Curie temperature is caused by Fermi-level-induced hole saturation. 
In order to provide a suitable treatment of the spin quasiparticle dynamics it is necessary 
to take into account the effects of disorder since the Mn ions are assumed to be distributed 
randomly with concentration c. This is positional disorder. There is variation of site-energy 
of nonmagnetic origin due to the substitution of A atom with Mn ion. The detailed nature 
of the disorder is not fully clear. In paper [36], it was shown that the dominant fraction 
of the Mn atoms are on either substitutional sites or specific sites shadowed by the host 
atoms. This reveals that the majority of the Mn atoms are on specific ( nonrandom ) sites 
commensurate with the lattice, but that does not necessarily imply that all of the Mn atoms 
are in substitutional positions. For x > 0.05 an increasing fraction of Mn spins do not 
participate in ferromagnetism. It can be related with an increase in the concentration of 
Mn interstitials accompanied by a reduction of T c . There are indications of an increase in 
Mn atoms in the form of random clusters not commensurate with the GaAs lattice. 
It follows from Eq.(94 ) that the spin dynamics of a modified KLM will be described by the 
GFs in the lattice site representation for a given configuration 



«S+\S j » «<7+|<7 j 



>> 



and instead of Eq.(28) the lattice GF should be considered 




(95) 
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In order to provide a simultaneous and self-consistent treatment of the quasi-particle dy- 
namics including the effects of disorder, a sophisticated description of disorder should be 
done. Most treatments remove disorder by making a virtual-crystal-like approximation in 
which the Mn ion distribution is replaced by a continuum. A more sophisticated approach 
for treating the positional disorder of the magnetic impurities inside the host semiconductor 
is the CPA [30]- [32]. The CPA replaces the initial Hamiltonian of disordered system by an 
effective one which is assumed to produce no further scattering [38] . It describes reasonably 
well the state of itinerant charge scattering in disordered substitutional alloys A\_ X B X . 
In order to simplify the discusson here, we will deal with a much simpler and less sophisti- 
cated description. The approximation discussed below should be considered as a first, crude 
approximation to a theory of disorder effects in DMS. Since the detailed nature of disorder 
in DMS is not yet established completely, we will confine ourselves to the simplest possible 
approximation. Let us remind that the IGF method is based on the suitable definition of the 
generalized mean fields [16]. To demonstrate the flexibility of the IGF method, we show 
below how the mean field should be redefined to include the disorder in an effective way. 
The previous definition of the irreducible spin operator, Eq.( 32), should be replaced by 

(Sg) ir = Sg — c< S z >5 q: o; (al +qa a pa ) ir = (ip +qa a pa — < a) pu a pa > 5 q ,o (96) 



Here < S z >= S z corresponds to the configuration average. The average < S z > denotes the 
mean value of S z for a given configuration of all the spins. We omitted here the variation 
of site energy of nonmagnetic origin. The consequences of this choice manifest themselves. 
It means precisely that in a random system the mean field is weaker as compared to a 
regular system. The approximation is conceptually as simple as an ordinary mean field 
approximation and corresponds to the virtual crystal approximation. The situation is then 
completely analogous to the previous one considered in the preceding sections. We get for 
the configurationally averaged GFs 

2cS? 



« St\S; »v =« S£\Sl k » ^ — r (97) 

*' 3 k 1 k uj - Im + 2I 2 cS z Xo(k,u) 
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Here S z = iV~ 1//2 < S z >. These simple results are fully tractable and are the reason for their 
derivation. 

It is worth to note that in the case of the modified Zener model which contains the correlation 
(Hubbard) term, the effects of disorder should be considered on the basis of a similar model 
[39] 

H = tiAaUjo + U J2 V i n i\ n H - J2 2Iu i°iSi (99) 

ij <y i i 

The Coulomb repulsion is assumed to exist only on lattice sites occupied at random by 
Mn atoms. The approach mostly used [39] to calculate stiffness constant within a random 
version of the Hubbard model was based on the random phase approximation, where the 
electron-electron approximation was taken into account in the Hartree-Fock approximation 
and the disorder in the CPA. It is therefore very probable that within this approach the 
formation of magnetic clusters can be reproduced; the formation of the clusters is thus 
strongly enviromental-dependent. However, the calculation of the spatial GF Eq.( 95), for 
the model, Eq.( 99), is rather a long and nontrivial task and we must avoid considering this 
problem here. We hope, nevertheless, that the description of the disorder effects, as given 
above, gives a good first approximation as far as the the irreducible Green functions method 
is concerned. A more detail consideration of the state of itinerant carriers in DMS including 
a more sophisticated treatment of disorder effects will be carried out separately. 

VII. CONCLUSIONS 

In summary, we have presented an analytical approach for treating the spin quasi-particle 
dynamics of the generalized spin-fermion model, which provides a basis for description of the 
physical properties of magnetic and diluted magnetic semiconductors. We have investigated 
the influence of the correlation and exchange effects for interacting systems of itinerant car- 
riers and localized spins using the ideas of quantum field theory for interacting electron and 
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spin systems on a lattice. The workable and self-consistent IGF approach to the decoupling 
problem for the equation-of-motion method for double-time temperature Green functions 
has been presented. The main achievement of this formulation was the derivation of the 
Dyson equation for double-time retarded Green functions instead of causal ones. That for- 
mulation permits to unify convenient analytical properties of retarded and advanced GF and 
the formal solution of the Dyson equation which, in spite of the required approximations 
for the self-energy, provides the correct functional structure of single-particle GF. The main 
advantage of the mathematical formalism is brought out by showing how elastic scatter- 
ing corrections (generalized mean fields) and inelastic scattering effects (damping and finite 
lifetimes) could be self-consistently incorporated in a general and compact manner. In this 
paper, we have confined ourselves to the elastic scattering corrections and have not consid- 
ered the damping effects. This approach gives a workable scheme for definition of relevant 
generalized mean fields written in terms of appropriate correlators. A comparative study of 
real many-body dynamics of the generalized spin-fermion model is important to characterize 
the true quasi-particle excitations and the role of magnetic correlations. It was shown that 
the magnetic dynamics of the generalized spin-fermion model can be understood in terms 
of combined dynamics of itinerant carriers, and of localized spins and magnetic correlations 
of various nature. These applications illustrate some of subtle details of the IGF approach 
and exhibit their physical significance in a representative form. 

As it is seen, this treatment has advantages in comparison with the standard methods of 
decoupling of higher order GFs within the equation-of-motion approach, namely, the follow- 
ing: 

At the mean-field level, the GF one obtains, is richer than that following from the stan- 
dard procedures. The generalized mean fields represent all elastic scattering renormalizations 
in a compact form. 

The approximations ( the decoupling ) are introduced at a later stage with respect to other 
methods, i.e., only into the rigorously obtained self-energy. 

The physical picture of elastic and inelastic scattering processes in the interacting many- 
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particle systems is clearly seen at every stage of calculations, which is not the case with the 
standard methods of decoupling. 

The main advantage of the whole method is the possibility of a self- consistent description 
of quasi-particle spectra and their damping in a unified and coherent fashion. 
Many results of the previous works are reproduced mathematically more simply. 
Thus, this new picture of an interacting spin-fermion system on a lattice is far richer and 
gives more possibilities for analysis of phenomena which can actually take place. In this 
sense, the approach we suggest produces more advanced physical picture of the quasi-particle 
many-body dynamics. Our main results reveal the fundamental importance of the adequate 
definition of generalized mean fields at finite temperatures, which results in a deeper insight 
into the nature of quasi-particle states of the correlated lattice fermions and spins. The key 
to an understanding of the situation in DMS lies in the right description of the interplay 
of interactions and disorder effects for coupled spin and charge subsystems. Consequently, 
it is crucial that the correct functional structure of generalized mean fields was calculated 
in a closed and compact form. The detailed consideration of the state of itinerant charge 
carriers in DMS along this line will be considered separately. 
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